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The electric dipole moments (EDMs) of heavy nuclei, such as ^^^Hg, ^^^Ra and ^^^Rn, can 
be enhanced by the Schiff moments induced by the presence of nearby parity-doublet states. 
Working within the framework of the maximally CP-violating and minimally flavour- 
violating (MCPMFV) version of the MSSM, we discuss the maximal values that such 
EDMs might attain, given the existing experimental constraints on the Thallium, neutron 
and Mercury EDMs. The maximal EDM values of the heavy nuclei are obtained with the 
help of a differential-geometrical approach proposed recently that enables the maxima of 
new CP-violating observables to be calculated exactly in the linear approximation. In the 
case of ^^^Ra, we find that its EDM may be as large as 6 to 50 x 10"^^ e ■ cm. 
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1 Introduction 



Electric dipole moments (EDMs) are among the most promising potential signatures for CP- 
violating physics beyond the Standard Model (SM), and one of the most promising options 
for extending the SM is provided by supersymmetry (SUSY) [1]. The minimal SUSY 
extension of the SM (the MSSM) already contains many possible CP-violating phases, 
even in its minimally flavour- violating (MFV) version. The maximally CP-violating MFV 
version of the MSSM, the MCPMFV model [21l3], has six CP-violating phases, to which may 
be added the QCD vacuum phase ^qcd- These phases are tightly constrained by the present 
experimental upper limits on the EDMs of ^°^T1, the neutron and ^^^Hg. Nevertheless, one 
could in principle envisage (accidental) cancellations |1] between the contributions to the 
measured EDMs of the six (seven) phases of the MCPMFV model (including 6'qcd), which 
might leave open the possibility of large net contributions to other EDMs. However, for 
any fixed values of the CP-conserving MCPMFV parameters, the compact ranges of the 
CP- violating parameters imply that the value of any other CP- violating observable, e.g., an 
EDM, is necessarily bounded. The question then arises whether the prospective sensitivity 
of any proposed experiment reaches below the maximum value attainable in any given 
theoretical framework, such as the MCPMFV model (with or without the possibility that 
^QCD 7^ 0). Clearly, any experiment that has insufficient sensitivity to search below the 
maximum value is not interesting for testing the MCPMFV model, whereas any experiment 
capable of reaching below the maximum value may either make a measurement or exclude 
part of the MCPMFV parameter space. 

In a recent paper [5] we proposed a novel analytical technique, based on a differential- 
geometrical construction (see also [B]), for finding the maximal values of CP- violating ob- 
servables subject to the existing EDM constraints, which is exact in the linear approxima- 
tion. We applied this technique to find maximal values of the EDMs of the Deuteron and 
muon, the CP- violating asymmetry in 6 — )■ S7 decay, Acp, and the Bg mixing phase [5]. We 
found that, whereas the EDM of the Deuteron in the MCPMFV model (allowing also for 
^QCD 7^ 0) might be one (two) orders of magnitude larger than the prospective experimental 
sensitivity, and Acp might also be detectable, the EDM of the muon and the contribution 
to the Bg mixing phase in the MCPMFV model are likely to be too small to be observable 
in the near future. 

In this paper we extend the applications of our analytical method [5l|6] to calculate the 
maximal values in the MCPMFV of the EDMs of some nuclei that are enhanced by the Schiff 
moment contributions associated with nearby parity-doublet states [7]. An experimental 
campaign is now being considered for HIE-ISOLDE that could search for EDMs of radium 
isotopes [S], accompanied by measurements of octupole collectivity in radium isotopes [H] 
that would be needed to interpret EDM measurements in terms of time-reversal violating 
interactions. 

We find that values of the ^"^^Ra EDM that are considerably larger than 10 e • cm 
(well within the estimated sensitivity of the proposed HIE-ISOLDE experiment) are possi- 
ble in the MCPMFV model. How much larger depends quite sensitively on the implemen- 
tation of the experimental constraint on the EDM of ^^^Hg. As discussed below, several 
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theoretical calculations of the EDM are available, and most of them give very similar 
allowed ranges for the EDM of ^^^Ra, in the range 6 to 10 x 10~^^ e ■ cm. However, one the- 
oretical calculation of the ^^^Hg EDM yields much weaker constraints on the CP- violating 
parameters of the MCPMFV, and hence yields a much larger maximal value of the ^^^Ra 
EDM, namely 50 x 10"^'' e ■ cm. 



2 SchifF Moments for Selected Nuclei 

A CP- violating atomic EDM may arise from intrinsic EDMs of the constituent nucleons and 
atomic electrons, CP-odd electron-nucleon interactions, and the CP-odd nuclear moment 
known as the Schiff moment. The Schiff moments of several atoms have been calculated in 
the literature, and expressed in terms of sums of products of CP-even and CP-odd pion- 
nucleon-nucleon [ttNN) couplings. We recall that the CP-even strong nNN interaction is 
given by [TO] 

^Tn'n = 9-nNN (pi-i^mT"^ + m75p7r") + {pi-f^p - ni-f^n) 7r° , (2.1) 
where q^nn = 13.45. On the other hand, the CP-odd (T-odd) nNN interactions are given 

by 



4nn = -qTmnNt'^Nti'^ + -g^^LNNrc' + (Nt'^Ntx^ - ?>Nr^NTx') 

= g^^NN (jj-riTT^ + npvr") + {pp - nn) 7r° 

+ 9nNN {pp + nn) 7r° 

\/2 iprni^ + npTX~] — 2 [pp — nn) 7r° (2.2) 



in terms of the isoscalar gfj^^-, isovector ^^tvat; ^'^^ isotensor ^^^^v T- violating pion-nucleon 
couplings. 

The Schiff moment is linear in the CP-odd nNN couplings g^jyfi^, and may be written 
as 112] 



S = [ao + b) g^^NNgnNN + (^1 gnNNgnNN + {(^2 — b) 

girNNgjiNN- 

(2.3) 

Here the coefficients specify the dependence of the Schiff moment on the CP-odd in- 
teractions and the coefficient b specifies its dependence on the nucleon dipole moments. 
The coefficients with i = 0,1,2 and b depend on the type of atom of interest, and some 
theoretical estimates for ^^^Hg, ^^^Ra and ^^^Rn, in units of e ■ fm are collected in Table [1] 

In the case of the ^^^Hg Schiff moment, the results [13] adopted in [1T1[TU[T5] are 
significantly different from the more recent calculations of [T2l[T6]. In particular, the coef- 
ficient Oo may be enhanced by a factor of ~ 30, whilst the coefficient ai may be reduced 
by a factor of ~ 10 and could even take the opposite sign. After considering the different 
numbers from the most recent calculations [12] in Table [H we have used the SIII (HF) 
calculation for our numerical illustration. As motivation, we note that the first SLy4 and 
the SV calculations do not yield the right quantum numbers for the ground state, that the 
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Table 1: The coefficients ai,i = 0,1,2 and b of the Schiff moments of^^^Hg, "^"^^Ra, and 
'^^^Rn, expressed in units of e ■ fm^. The labels HB and HFB stand for calculations in the 
Hartree-Fock and Hartree-Fock-Bogoliubov approximations, respectively: see Ref. fW^ for 
details. We have changed the signs of the coefficients oq, ai, and b to follow the conventions 
of Ref m^. 



Atom 



Ref. Interaction 



-ao 



-ai 



a2 



-b 



199 



Hg 



0.0004 



0.055 0.009 



SkO' 

(average) 



0.010 
0.007 



0.074 
0.071 



0.018 
0.018 



SLy4 (HF) 
SIII (HF) 
SV (HF) 
SLy4 (HFB) 
SkM* (HFB) 



0.013 
0.012 
0.009 
0.013 
0.041 



-0.006 
0.005 
-0.0001 
-0.006 
-0.027 



0.022 
0.016 
0.016 
0.024 
0.069 



0.003 
0.004 
0.002 
0.007 
0.013 



225 



Ra 



[T7] SkO' (zero range) 



-5.1 



10.4 -10.1 



PJ SkO' 



-1.5 



6.0 



-4.0 



211 



Rn 



il9j 



0.019 



-0.061 0.053 



12] SLy4 
SkM* 
SHI 



0.042 
0.042 
0.034 



-0.018 
-0.028 
-0.0004 



0.071 
0.078 
0.064 



0.016 
0.015 
0.015 



second SLy4 calculation was found to be unable to project HFB states, and that the SkM* 
calculation gives a spherical minimum. 

The CP-odd couplings g^jfi^ may be generated by chromoelectric dipole moments 
(CEDMs) of the quarks and/or dimension-six four-fermion interactions. In Ref. [20], one 
may find full expressions for the contributions to the CP-odd isoscalar g^)fj<^ and isovector 
g^^iy couplings from the CEDMs of the light quarks E|: 

g^:L = (-0.5 to 1.5) X 10--^^-"^ + ^^)/^^ l^^^'^^l 



10-26 cm (225MeV)3' 

W - ^-iXlO 10-26 cm (225MeV)3- ^^'^^ 

We note that a 'best value' is not available for g^^l^j^- For definiteness, we have taken 
9wNn/9^nn = 0-2 {du + dd)/{du — d^), as follows from assuming that the couplings g^^N 
and g^NN proportional to the matrix elements of uu — dd and uu + dd in the nucleon 
state, respectively. It is known that the matrix element with the minus (— ) sign is smaller 
than the one with the plus (+) sign by a a factor ~ 5 to 10 [21] . 



'We adopt the conventions and notations of [MjlTS 
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We have also included the contribution of the dimension-six four-fermion interactions 
to the isovector coupling which can be enhanced for a large value of tan (3 



9^L = ^{N\Cdddd + CsdSS + CMbb\N) 



-8 X lO^^GeV^ 



h 3.3k hi — 0.25k 



where the couplings appearing in (12. 5 p are defined via the interaction Lagrangian 



(2.5) 



Af = EC//'(/7)(/'^75/'), (2.6) 
/,/' 

and K = {NlrrisSslN) /220 MeV ^ 0.50 ± 0.25. The contribution of the four-fermion 
interactions to the isoscalar coupling g^^j^^ is ignored and the isotensor coupling g^p^N, 
which changes the isospin by two units, is neglected in this work (as in Ref. |20]), since it 
can be generated only at the expense of an additional rriu — rrid suppression. 

In the following, we concentrate on the EDMs of ^^^Hg and ^^^Ra atoms, because 
there is already a stringent experimental upper limit on the former, and there is a pro- 
posal to measure the latter at HIE-ISOLDE [8]. Moreover, we note that the estimates of 
the contributions to the Schiff moment of ^^^Rn shown in Table [1] are considerably less 
favourable. 

The experimental constraint on the EDM of ^^^Hg [211I2S], 

dij^l < 3.1 X 10-29e-cm (95%C.L.), (2.7) 

is among the most stringent constraints on CP violation imposed by atomic EDM experi- 
ments. There are several calculations of the ^^^Hg EDM in terms of the Schiff moment in 
the literature, which can be cast in the general form 

du,[S] = 10-17 C|g e ■ cm X (^-^-g) , (2.8) 

with the coefficient taking the values C^g = -4 [2H1[27], -2.8 |2B] and +5.07 [22]. We note 
that we have changed the sign of the most recent calculation based on relativistic coupled- 
cluster theory [29], so as to match the conventions in the previous calculations [26l - [28] 
based on coupled perturbed Hartree-Fock calculations, and that it has a different sign. For 
our numerical calcluatons we take the estimates of the coefficient C^g that have negative 
signs. To be specific, we have used the following four calculations of the Mercury EDM: 



From Ref. [TT 

• Taking = — 2.8 [28] and the coefficients from Ref. [13 



dUS] ~ 1.8 X 10-3 e^S^/GeV, (2.9) 



d''[S] ~ 7.6 X 10-^e^S^/GeV + 1.0 x 10-^ e^S^ /GeV , (2.10) 
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• Taking C^^ = -2.8 [2H] and the average coeffcients from Ref. [T6] . 

diliS] ~ 1.3xlO-^e^S^/GeV + 1.4xlO"=^e^S^/GeV, (2.11) 

• Taking C^^ = -2.8 [28] the SIII(HF) coeffcients from Ref. [12], 

d^^[S] - 3.1 X lO-^e^S^/GeV + 9.5 x 10"^ e /GeV . (2.12) 

Compared to the expression ^Hgl'^'] used in Refs. [HllIS], d^g[S] is shghtly smaller, and has 
a negligible isoscalar contribution, which is actually ignored in ci^gl'S']. Whilst the isoscalar 
contribution is still small in c?^g[5'], it can be larger than the isovector contribution by a 
factor ~ 3 in d^g[S]. However, we observe that I'^Hgl'S']! is smaller than |(iHg['S']| by about 
an order of magnitude. 

On the other hand, the EDM of ^^^Ra is related to its Schiff moment by [28] 

dR.[S] = -8.5x 10-^'e-cmx (^— ^) . (2.13) 

Taking the coefficients given in Ref. [IB] , see Table [H we obtain 

dn^S] - -8.7x 10-2 e^S^/GeV + 3.5 X 10^1 e^Sjv/GeV. (2.14) 

We note that the g^^lu^ contribution the EDM of ^^^Ra is about 200 times larger than to 
the Mercury EDM rfHgfS"], an enhancement due to the existence of a nearby parity-doublet 
states [T7] . 



3 DifFerential-Geometrical Optimization Method 

We briefly review our powerful analytical approach for flnding the optimal choice of CP- 
odd phases which maximize the size of a given CP- violating observable O, while remaining 
compatible with the present EDM constraints [6] . We have applied this approach previously 
to estimate maximal values of the Deuteron and muon EDMs, the CP asymmetry in 6 ^ 57, 
and the phase in Bg mixing [5], and it may be applied similarly to the case where the 
observable O is the EDM of ^^^Ra. 

We consider a theory such as the MCPMFV SUSY model with six CP-odd phases, 
represented by a 6D phase vector, subject to three EDM constraints denoted by E"'''^''^ = 0, 
corresponding to the non-observation of the Thallium, neutron and Mercury EDMs. For 
any given value of the CP-conserving parameters in the MCPMFV model, we may expand 
these EDMs and the observable O in the small-phase approximation, deflning the four 6D 
vectors E"'*''^ = V-E"'*''^ and O = VO, and we assume that the four vectors E"''''^ and O 
are linearly independent. 

We then introduce the triple exterior product 



6 



where the Greek indices label the components of the vectors in the 6D space, i.e., a, /3, 7 = 
1,2, ...,6. The square brackets on the RHS of f l3.15p indicate that the tensor Aa/s-y is 
obtained by fully antisymmetrizing the vectors E^, and E^ in the indices a,/3,7, i.e., 
A-ajS-y = —A.i3a^ = —Aa^j^, etc. Borrowiug a term from the calculus of differential forms, 
Aaji^ is a 3-form. We also introduce the 2-form 



OxE;E'^E% (3.16) 

where summation over repeated indices is implied and e^^Xpar is the usual Levi-Civita ten- 
sor generalized to 6D. In the language of differential forms, Bf^^, is, up to an irrelevant overall 
factor, the Hodge-dual product between the 1-form Ox, representing the CP-violating ob- 
servable, and the 3-form Aap-y- 

The components $* of the optimal EDM-free direction maximizing O can now be 
obtained from the Hodge-dual product of the 3-form Aj^^s and the 2-form B^y. Explicitly, 

^al3"i5pu Ap^fS — M Eal^-ySpu^puXpar E'j!j E^^ E^ OxEp E^ E^ , (3-17) 

where we have included an unknown overall normalization factor A/". By construction, the 
6D phase vector $* is orthogonal to the three vectors E'^'''''^, and therefore satisfies the 
desired EDM constraints, E"" = E'^ = E'^ = 0, in the small-phase approximation. We 
observe that the magnitude 0* = |$*|, and hence the overall normalization factor A/", can 
only be determined by a numerical analysis of the actual experimental limits on the three 
EDMs. As in the 3D example, the maximum allowed value of the CP-violating observable 
O is given in the small-phase approximation by 

O = 0* $: = ±Ar \e^ua^^s e^^xpar Ox E^ E'^ E^ E; E^ E^, \ , (3.18) 

where the caret denotes the components of a unit-norm vector. As discussed in [5], 
quadratic and higher-order derivative terms with respect to the CP-odd phases will gener- 
ically prefer a particular sign for the optimal value of O. 

We can also allow for the possible presence of a non-zero strong CP phase ^qcd in the 
theory, in which case the corresponding CP-odd phase vector $ becomes seven-dimensional 
(7D) in the MCPMFV SUSY model. The generalization of the above construction of the 
optimal value of the observable O is discussed in Section 5 of [5]. 



4 Results 

In this Section we use the above the differential-geometrical approach to analyze the maxi- 
mal value of the ^^^Ra EDM obtainable in CP- violating variants of the following represen- 
tative CMSSM benchmark scenario which predicts the mass spectrum of SUSY particles 
in the sub-TeV region: 

|Mi,2,3| = 350 GeV, 

Ml^ = Ml^ =Ml = Ml = Ml = Ml = Ml = (100 GeV)^, 

\A^\ = \Ad\ = \Ae\ = 100 GeV, (4.19) 
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at the GUT scale, introducing non-zero CP- violating phases and varying tan/3 (Msusy)- 
We adopt the convention that = 0°, and we vary independently the following six 
MCPMFV phases at the GUT scale: $i, $2, $3, $Au) "^a^, and ^a^- In the 7D case, in 
addition to the 6 CP phases, we consider a non-zero strong CP phase 6'qcd- The scenario 
fl4.19p considered in this work is similar to one with | Mi 2,3] = 250 GeV that we considered 
previously |2l[5l[14]. The somewhat larger value |Mi^2,3| = 350 GeV is chosen here for 
consistency with the recent SUSY search results reported by the CMS Collaboration at the 
LHC [30]. When tan/3 = 10, $1,2,3 = 0°, and ^A^AdAc = 180°; the previous scenario with 
1^1,2,3] = 250 GeV became the well-known SPSla point [31] (also known as benchmark B). 
The scenario f l4.19p is more similar to benchmarks C, G and I of [31] when tan/3 = 10, 20 
or 35, respectively. Our calculations of the EDMs are based on Refs. [5l[T3], which include 
the two-loop diagrams mediated by the 'j-H'^-W^ and ■y-W'^-W'^ couplings, and we take 
into account the effects of the different computations of the Schiff moment of the Mercury 
nucleus as explained in Section [2l 



4.1 The MCPMFV SUSY model with 6 CP phases 

In order to analyze the scenario f l4.19p . we first make Taylor expansions of the following 
four EDMs in terms of the MCPMFV CP-violating phases: 

dTi/d'ir , rfn/d^^P , dnjdlf , dnj {10~'' e ■ cm) , (4.20) 

where we choose the following normalization factors: 

4fP = 9xl0"2^e-cm, df^^ = 3x lO"'*' e-cm , d^f = 3.1 x 10"'^ e-cm , (4.21) 

which correspond to the current experimental limits on the EDMs of Thallium [32], the 
neutron [33], and Mercury [2^ [25]. The normalization factor for the ^^^Ra EDM, namely 
10~^^ e ■ cm, is typical of the estimated experimental sensitivity. 

In Fig. [H we show the absolute values of the components of the three 6D MCPMFV 
vectors characterizing the existing EDM constraints and the 6D vector representing the 
dua observable, for the scenario (14.19^ varying tan /3 in a small-phase expansion around 
the CP-conserving point $1 = $2 = '^'s = ^a^ = = ^a^ = 0° Q. The solid lines 
represent the CP-violating phases of the gaugino mass parameters, and the dashed lines 
the trilinear A parameters, respectively. The components corresponding to the CP- violating 
gaugino phases dominate in all cases, increasingly as tan /3 grows, with the exception of 
the $3 component of E'^'^'. For example, when tan /3 = 40, \dT\\ and \dn\ are larger than 
the current experimental limits by factors ~30 and ~5, respectively, even when $2 = 1°, 
whereas \d^g\ is larger than the current limit by a factor ~15 (~8) when $3 = 1° ($2 = 1°) 
0. We see that |(iRa| can be as large as ~ 100 x 10~^''ecm (~ 40 x 10~^^ecm) if $3 = 1° 
($2 = 1°). 

^We also display (as magenta lines) the corresponding components for 7th components corresponding 
to the QCD phase 9, which we discuss later. 

Throughout this work, the CP- violating phases are measured in units of degrees. 
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10 20 30 40 10 20 30 40 

tan[3 tanP 



Figure 1: The absolute values of the components of the three vectors Eq, of MCPMFV CP- 
violating phases representing the present EDM constraints on d^x (upper left), (upper 
right), and c/ng (lower left), and those of the vector representing the EDM o/^^^Ra (lower 
right) in small-phase expansions around the CP-conserving point, as functions of tan P for 
the scenario (^[7^. We use the computation employing the QCD sum rule technique and 
the estimate d^^ h2.9) for the neutron and Mercury EDMs, respectively. The CP-violating 
phases are measured in degrees. 

The relative contributions of the different CP- violating MCPMFV phases to (ing vary 
according to the choice of theoretical calculation, as shown in Fig. Specifically we 
observe the $3 contribution exhibits strong variations, and that it is much suppressed if 
the I^Hgl calculation is used. 

In Fig. |3]we show the cosines of the angles between the 6D MCPMFV vector repre- 

^We again display (as magenta lines) the corresponding components for possible 7th components cor- 
responding to the QCD phase 6, which we discuss later. 
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10 20 30 40 10 20 30 40 

tanP tanP 

Figure 2: The absolute values of the components of the 6D vector of MCPMFV CP- 
violating phases representing the present EDM constraints on the Mercury EDM, as ob- 
tained using the four different calculations of d^g discussed in Section\^ d^^ (upper left), 
d^g (upper right), d^^ (lower right), and d^^ (lower left). The line styles are the same as 
in Fig. [II The CP-violating phases are measured in degrees. 



senting the observable 0'^^=' and the EDM vectors, which are defined by 

^o,E^ = i^pI^ ■ (4-22) 

These quantities display the degrees of degeneracy, or alignment, between the observable 
0"'R- and the EDM vectors E". We see that the alignments between d^^ and d^i and d^^ 
are quite small. On the other hand, the alignments between c^Ra and c^Hg^'^^^ larger, 
though decreasing as tan /3 grows. This can be understood by observing that the additional 
contributions to the Mercury EDM, beyond those from the couplings g^pl:J'\ increase with 
tan (3. 
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Figure 3: The cosines of the angles between the observable 6D vector O and the EDM- 
constraint 6D vectors E'^ti fi)iack solid), E'^" (red solid), E'^hs (blue solid), E'^Hg 
dashed), E (blue dotted), and E (blue dash-dotted) as functions o/tan/3. 

Having the vectors representing the EDM constraints and the observable rf^a in hand, 
we now combine them to construct the optimal directions in the 6D space of CP-violating 
MCPMFV phases, using f l3.17p . so as to maximize d-^a, in the linear approximation. For 
comparison, we also consider two reference directions, which have $i = $Ae = and 
$2 = '^'s = 0, respectively. These two reference directions can be constructed by defining 

K = ^feap,s,u El E\ El iV(^) iV(2) ^ (4_23) 

where, for each direction, the two null directions N^'"^^ are chosen as 

iV« = (1,0,0,0,0,0), iVf = (0,0,0,0,0,1) (4.24) 

for the direction $i = = 0, and 

TVW = (0,1,0,0,0,0), TVf ^ (0,0,1,0,0,0) (4.25) 
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A* A* A* A* A* A* 




Figure 4: The absolute values of the six components of the normalized vectors in the di- 
rections where the size of dua. is optimized, as obtained using the calculations d^^ (upper 
left), d^g (upper right), d^^ (lower right), and d^^ (lower left) for the Mercury EDM. The 
CP-violating phases are measured in degrees. 



for the direction ^2 = ^3 = 0. 

We display in Fig. S] the absolute values of the six components of the normalized 
optimal vectors in the direction along which d^^a, is maximized, as obtained using d^^ (upper 
left), rf^g (upper right), d^^ (lower right), and d^^ (lower left) for the Mercury EDM. We 
first observe that the $1,2,3 components (solid lines) are relatively small, and decrease as 
tan /3 increases. Hence, all the optimal directions are mostly given by some combination 
of (black dashed line) and (i'gcI dashed line) directions implying, for tan /3 = 40, 
that ($A„,Ad)™^'' ~ 0* whereas ($1,2,3)™'^'' ~ 0* x 10"^, as will be shown in the following. 

In Fig. we consider the products $* ■ O of the 6D vectors in the normalized optimal 
directions for dua, and the Radium EDM, taking account of the uncertainty of the Mercury 
EDM calculation. The products determine the sizes of d^a, along its optimal direction 
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Figure 5: The products $* ■ O (black solid) for the optimal d-Ra direction, as obtained using 
the calculations d^^^ (upper left), d^^ (upper right), d^g (lower right), and d^^ (lower left) 
for the Mercury EDM. For comparison purposes, we also show the products along the two 
reference directions with = (red dashed) and $2,3 = (blue dotted). 



through the relations given in (13.181) when (p* = 1°. As shown below, (p* could be as 
large as ~ 100° before the small-phase approximation breaks down and one of the three 
EDM constraints is violated. We observe that the direction constructed using the geometric 
prescription given in Section [3] indeed gives the larger values of ciRa than do the two reference 
directions with $i,yi^ = and $2,3 = 0. We note in particular that, in the case of d^g, 
the Radium EDM may become about an order of magnitude larger than is possible with 
another choices for the Mercury EDM calculation. 

Fig. |6] shows the Thallium (solid), neutron (dashed) and Mercury (dotted) EDMs 
along the directions chosen to optimize c^Ra for the four choices of dug when tan /3 = 40. 
We see that (p* could be as large as about 100°, with the most important constraints being 
provided by the Thallium (Mercury) EDM in the cases where the c^Hg^^'^^ ('^Hg) used. 
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Figure 6: The values of the three EDMs along the directions optimized for c/Ra? obtained 
using d^^ (upper left), d^^ (upper right), d^^ (lower right), and d^^ (lower left) for the 
Mercury EDM. The ratios dTi/d^f^ , d^/di^^ , and d-a^/d^^ are denoted by the black solid, 
red dashed, and blue dotted lines, respectively. The scenario ^4-^9 ) is assumed, with the 
choice tan/3 = 40. 



respectively. 

In Fig. [TJ we show the maximal values of c^Ra attainable in the 6D case after imposing 
the three EDM constraints. We see the that small-phase approximations for (i^g^'^^^ break 
down for 0* ~ 40°, as seen by comparing with the calculations along the two reference 
directions, and that d-^^. can be as large as ~ 6 x 10~^^ e ■ cm or more. On the other hand. 
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e ■ cm. 



in the case of d^^, rf^a can be as large as ~ 50 x 10" 

Finally, in Fig. |H] we show the 6 CP-violating phases at the GUT scale (top and 
middle) and the 3 CP-violating phases of the third-generation A parameters at the SUSY 
scale (bottom). We observe that the CP phases of the gaugino mass parameters $i, $2 and 
$3 can only be as large as 2°, 0.4° and 0.4°, respectively, whereas (^a^) at the GUT 
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Figure 7: The EDM of '^^^Ra in units of 10"^"^ e cm for the scenario Ili4-19\ ) with tan /3 = 40 
along the directions optimized for d^^. (black solid) taking d^^ (upper left), d^^ (upper right) , 
^Hg (lower right), and d^^ (lower left) for the Mercury EDM. We have imposed the EDM 
constraints |(iTi/c?ff^l — 1? \dn/d^^^\ < 1, and Mng/c^Hg^l — 1- ^or comparison, we also 
show djia along two reference directions with the choices $i,Ae = (''"^d dashed) and $2,3 = 
(blue dotted). 



scale can be as large as ~ 100° (60°), as previously seen in Fig. HI These CP- violating phases 
are suppressed at the SUSY scale by RG running from the GUT scale [2], but sizeable 
non-trivial CP- violating phases are still allowed at the SUSY scale: 1$^^^ — 180° | ~ 2°, 
1$^, - 180°| ~ 4°, and - 180°| ~ 30°. 
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Figure 8: The 6 CP-violating phases at the GUT scale (upper and middle) and the 3 CP- 
violating phases of the third- generation A parameters at the SUSY scale (lower) along the 
directions optimized for d-Ra, obtained using d^^ (solid black), d^^ (red dashed), d^^ 
(blue dotted), and d^^ (magenta dash-dootted) for the Mercury EDM. We have imposed the 
EDM constraints |c?Ti/ciff^| < 1, K/rf^^^l < 1, and \dng/d^Y\ ^ 1- scenario U^^ j 
is assumed, with the choice tan/3 = 40. 



4.2 The 7D Case of non-zero 6'qcd 

Hitherto, we have imphcitly assumed that the CP-violating QCD ^-term: 

C 



cy — ~ 



(4.26) 



is neghgible, where G^'^'"' = e^'^P"G'^^^/2 and the parameter 9 is given by the sum of the 



QCD ^QCD and the strong chiral phase for the quark mass matrix: 

e = 0QCD + ArgDetMg. 



(4.27) 
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In the weak basis where ArgDet Mq = 0, we have 6 = Oqqd- 

The dimension-four operator f l4.26p would in general contribute to the neutron, Mer- 
cury and Radium EDMs, e.g., through the CP-odd pion-nucleon-nucleon interactions (12. 2p . 
Explicitly, for the neutron EDM, we use the estimate [5|[TT] 

dn{e) ~ 2.5 X 10"^^^ e- cm. (4.28) 

For the ^-induced Mercury EDM, we neglect the contribution from the g^)fj\f coupling and 
use 

4f '"^'^^(^1 - iC^r''''' X 10-^ GeV'^)etU0) , (4-29) 
where C^^ = 1.8, = 1.0, C^^^ = 1.4 and C^^ = 9.5 x 10"^, with 

£lj,{e) ^ 1.1 X 10-3^. (4.30) 

Finally, for the Radium EDM, we use 

dU^) ~ (3.5 X 10-1 GeV-')e£l^{e) . (4.31) 

Henceforward, we normalize 6 in units of lO^^^: 

e = 9x 10^°. (4.32) 

With this normalization, when 6' = 1, we have dn{0) = 2.5x10"^^ e-cm, which is very close to 
the current experimental bound d^^^ = 3 x 10"^^ e-cm, and we find (i^a — 0.8 x 10"^'' e-cm. 

We have analyzed the possible maximal values of in this 7D case including 6 fol- 
lowing a procedure similar to that we used in the 6D case. Looking again at Fig. [H we note 
the horizontal magenta lines representing the 7th components of the vectors representing 
the present EDM constraints on o?ti (upper left), dj^ (upper right), and rfng (lower left), 
and of the vector representing the EDM of ^^^Ra (lower right). We observe that the 7th 
component is missing in the rfxi case, because this observable has no contribution from the 
6 term in our approach. On the other hand, the 6 component is close to unity for d^ (cf. the 
discussion at the end of the previous paragraph), ~ 10"^ for d^^, and somewhat less than 
unity for d^a- This implies that a measurement of d^^^ at the level of 10"^^ e cm would 
already be a competitive measurement of 6, even in the absence of the other MCPMFV 
phases. 

As seen in Fig. [21 we have made similar analyses using the d^g^^i'^^ calculations, 
finding similar 6 components in the first two cases, but a value about an order of magnitude 
smaller in the d^^ case. We have also analyzed (not shown) the cosines of the angles between 

,1,11, III, IV 

the observable vector O^^^ and the EDM-constraint vectors E E " and E Hg in 
the 7D model, as functions of tan/3. We find results that are very similar to the 6D case 
shown in Fig. 121 the most significant difference being quite small and limited to tan /3 < 10 
in the d^ case. 

We turn now to Fig. [9l which is the analogue of Fig. |H but including the QCD 6 term 
(j4.26p . We see significant differences at intermediate and large tan/3, where we note that 
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Figure 9: ^45 in Fig. but including the QCD 9 term ^.20^ . The 7th components of the 
normalized 7D optimal vectors are denoted by the solid magenta lines. 

the optimal vectors in all four models for d^g exhibit not only relatively large 6 components, 
but also larger components for the gaugino mass phases than in the 6D case. We have also 
analyzed (not shown) the products $* ■ O for the optimal c^Ra direction in the 7D space, 
as obtained using the calculations c/^g^'^^^'^^ for the Mercury EDM, finding results that are 
generally very similar to those shown in Fig. |5]with a little rise around tan/3 = 40. 

Figure [10] displays the values of the Thallium, neutron and Mercury EDMs along 
the directions optimized for c^Ra, analogously to Fig. [6] but this time in 7D. As before, the 
different panels display results obtained using d^^ (upper left), d^^ (upper right), (i^g (lower 
right), and d^ (lower left) for the Mercury EDM, and the ratios d-n/d^f^, d^/d^ , and 
dug/d^^ are denoted by the black solid, red dashed, and blue dotted hues, respectively. 
Also as before, the scenario (14.191) is assumed, with the choice tan /3 = 40. Comparing with 
Fig. O we see that the magnitudes of the Thallium and Mercury EDMs increase quickly as 
(j)* deviates from and it could be as large as only about 50°. 
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Figure 10: As in Fig. O but including the QCD 6 term U-^Gj )- 



Turning to the resulting 7D estimates of c^Ra shown in Fig. [TTl we see that they 
have somewhat smaller magnitudes than the 6D estimates shown in Fig. [7] in the linear 
approximation. However, larger values of the CP-violating gaugino phases are allowed in 
the 7D case than in the 6D case, as seen by comparing Figs. [S] and [T21 Specifically, we 
observe that $i, $2 and $3 could be as large as ~ 15°, ~ 4° and ~ 1°, respectively. We 
also see in Fig. [13] that sizeable 6 could be much larger than the upper limit of ~ 10^^*^ 
usually quoted, with values as large as 6' ~ 2.5 x 10"^ becoming possible in the presence of 
non-zero MCPMFV phases. 



5 Conclusions 



In this paper we have extended our previous analyses of the MCPMFV model, with its 
6 CP-violating phases to determine (in the linear approximation) the largest value of d-Ra, 
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Figure 11: As in Fig. but including the QCD 6 term U-^Gj )- 



that is allowed by the present constraints on the neutron, Thallium and Mercury EDMs, 
using the differential-geometric approach developed in [5]. Numerically, we obtain rather 
similar results whether we include the CP-violating QCD vacuum phase 9 in the analysis, 
or not. 

The results are much more sensitive to the theoretical treatment of the Mercury 
EDM constraint, and we compare the results obtained with four different calculations of 
dug. Three of them yield quite similar results for dua, but one calculations indicates a 
smaller dependence of d^g on the CP-violating phases of the MCPMFV model, and hence 
allows larger numerical values of these phase and, in general, larger values of dpta become 
possible. The maximal values we find for dpta using three of the d-^g calculations are typically 
~ 6 X 10~^^ e ■ cm or more, whereas the fourth calculation allows d^g, ~ 50 x 10~^^ e ■ cm. 

For comparison, we recall that there is a proposal to measure c/Ra with an sensitivity 
approaching ~ 10"^'' e ■ cm in one day of data-taking. This experiment would clearly have 
interesting potential to probe regions of the MCPMFV parameter space that have not been 



20 




explored by previous EDM experiments. 

This potential surely extends to many other models with several sources of CP viola- 
tion, which could also be analyzed using the differential-geometric approach [S] exploited 
here. As long as a limited number n of EDMs have been bounded (or measured) by ex- 
periments, any model with N > n CP-violating parameters will be underconstrained, and 
(partial) cancellations [1] are possible that would allow large values for other CP-violating 
observables. The MCPMFV model is one such example, in which = 6 (or 7 if 6' is in- 
cluded) and 77, = 3, so far. As we have shown in this paper using this differential-geometric 
approach, constraining or measuring c^Ra at the level of ~ 10"^^ e ■ cm or better would be 
a valuable addition to the existing arsenal of experimental probes of CP violation. 

However, it would still not complete the set of constraints needed for a model with 
A^ > 4, such as the MCPMFV. For this reason, other measurements, e.g., of the CP- 
violating asymmetry in 6 —j- 57 decay. We have discussed elsewhere the maximal value 
that this observable might take in the MCPMFV, and it would be an interesting and 
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Figure 13: The values of 6 = 6 x 10^° along the directions that maximizer in the linear 
approximation, as found by using the calculations d^^ (solid black), d^^ (red dashed), d^g 
(blue dotted), and d^^ (magenta dash-dootted) for the Mercury EDM. We have imposed the 
EDM constraints \dTi/d^f^\ < 1, \dn/df^^\ < 1, and {dng/d^Yl ^ 1- scenario gJ^ 
is assumed, fixing tan /3 = 40. 



complementary measurement to that of dua- 

We conclude by drawing reader's attention to potential caveats in the physics of CP 
violation. The baryon asymmetry of the Universe is one of the strongest pieces of evidence 
for physics beyond the Standard Model, as it cannot be generated successfully within 
the standard Kobayashi-Maskawa model of CP violation. There must be new sources of 
CP violation beyond the Kobayashi-Maskawa phase, and it behooves experiments to chase 
down all those within reach. Some CP- violating phases may manifest themselves at the TeV 
scale and be accessible to contemporary collider experiments, e.g. at the LHC. However, 
baryogenesis could equally well be achieved via CP-violating phases appearing at higher 
energy scales, and EDMs have the potential to probe beyond the TeV scale, in particular 
because the Standard Model Kobayashi-Maskawa predictions for EDMs are quite small. As 
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our analysis exemplifies, new EDM observables probe complementary region of parameter 
space including the strong CP phase ^qcd- Therefore, constraining or even measuring d^is. 
is an interesting experimental objective. 
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